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We report the possibility of completely destructive interference of three indistinguishable photons 
on a three port device providing a generalisation of the well known Hong-Ou-Mandel interference 
of two indistinguishable photons on a two port device. Our analysis is based on the underlying 
mathematical framework of SU(3) transformations rather than SU(2) transformations. We show 
the completely destructive three photon interference for a large range of parameters of the three 
port device and point out the physical origin of such interference in terms of the contributions from 
different quantum paths. As each output port can deliver zero to three photons the device generates 
higher dimensional entanglement. In particular, different forms of entangled states of qudits can be 
generated depending on the device parameters. Our system is different from a symmetric three port 
beam splitter which does not exhibit a three photon Hong-Ou-Mandel interference. 


INTRODUCTION 

The Hong-Ou-Mandel (HOM) effect (TJ [2], i.e., the 
completely destructive interference of two independent 
but indistinguishable photons, brought a paradigm shift 
to the field of quantum optics. Until the demonstration of 
the HOM effect the interference of independent photons 
was considered to be impossible. Such an interference 
effect manifests itself in the study of photon correlations 
rather than in intensity measurements. More specifically, 
if two single photons are sent from two different ports of 
a 50/50 beam splitter then the number of coincidence 
events at the two output ports vanishes. This follows 
from the fact that if two photons are indistinguishable 
with respect to their wavelength and polarisation and 
their wave packages overlap in time then the two differ¬ 
ent quantum paths interfere so that the two photons will 
never leave the beam splitter at different ports. If one 
of these parameters is changed the photons become dis¬ 
tinguishable and the dip in the observed coincidence rate 
starts to disappear. 

The effect is quite versatile and has been observed in a 
very wide class of systems. Besides systems with discrete 
optical elements like beam splitters it has been studied 
in other optical setups such as in integrated devices like 
evanescently coupled waveguides [SU and coupled plas- 
monic systems mi- It has also been studied in the ra¬ 
diation from two trapped ions uni , atoms HEHEI] , quan¬ 
tum dots m ii9 and for two different kinds of sources 

nsnsi- 

Since the original work of HOM one also has examined 
the kind of interference that can take place if two single 
photons are replaced by say two photons on each port or 
say by one at one port and two at the other port. Hereby 
one has found very interesting quantum interference ef¬ 
fects depending on the beam splitter reflectivity fTT)T(2T] . 
Another interesting possibility occurs if n photons arrive 


at each port of a 50/50 beam splitter - in this case the 
output ports never have odd numbers of photons [22]. 

In this letter we report a three photon interference ef¬ 
fect which is in the original spirit of the HOM effect - we 
examine the completely destructive interference of three 
indistinguishable photons on a three port device. We 
thus shift the focus from a two port device to a three 
port device. This brings a key change to the under¬ 
lying mathematical framework as we work with SU(3) 
transformations rather than SU(2) transformations. We 
specifically examine a three port integrated device con¬ 
sisting of a small array of three single mode evanescently 
coupled waveguides as these are relatively easy to fabri¬ 
cate [23[ |2l| ■ Although we tailor our discussion to cou¬ 
pled waveguide systems the results will be applicable to 
a wide class of bosonic systems described by the Hamil¬ 
tonian ([l]). For the three port device we have found an 
analytical expression for the completely destructive three 
photon interference. Thereby we produce a variety of two 
and three photon entangled states at the output ports. 

Our three port network is different from the symmet¬ 
ric multiport beam splitter which has been extensively 
studied for HOM like interferences [ 25 ] [ 25 ], as for the 
three port system such a splitter does not exhibit a per¬ 
fect three photon HOM interference. On the other hand, 
Campos m using the idea of Reck et al. [ 25 ] con¬ 
structed a SU(3) transformation involving beam splitters 
and phase shifters which leads to three photon HOM 
interference. Further, Tan et al. [29] showed how the 
SU(3) transformation involving beam splitters and phase 
shifters can lead to perfect photon interference depending 
on specific values of the parameters of SU(2) transforma¬ 
tions. While the papers of Campos m and Tan et al. 
[ 25 ] concentrate on configurations using beam splitters 
we investigate integrated devices. Such integrated de¬ 
vices are now used by many experimental groups. Studies 
of multi photon interferences in integrated devices have 
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the waveguide mode’s frequency is matched to the fre¬ 
quency of the input field. The coupling strength is given 
by the coupling coefficients g\ (between the first and sec¬ 
ond mode) and g 2 (between the second and third mode) 
and is basically determined by the distance between the 
guides. The calculation of the coupling parameters is text 
book material and can be found, e.g., in m■ Note that 
the coupling between the two outer modes is negligible 
for our geometry. 

The interaction Hamiltonian for this system reads 


FIG. 1. Scheme of a 3 x 3 waveguide array with continuous 
coupling between the inner mode and the outer modes 


been reported, e.g., in 0013SM]- Note that the ex¬ 
amination of multiport systems is also important in the 
context of Bosonic sampling [33ti36] [38| where the output 
distribution of Fock input states is sampled what is ex¬ 
ponentially hard to predict with classical computations. 

While our investigated setup is similar to the experi¬ 
ment studied in [39] . which relies on a 3D geometry in 
order to couple all three modes to each other, we will use 
a more simple 2D structure, where the outer modes are 
coupled to the inner mode but not to each other. With 
the setup of [35] it is possible to suppress [25] states of the 
form 12,1,0), which contain two, one and zero photons 
in the different output modes, but the output state will 
still contain a |1,1, l)-term corresponding to the coinci¬ 
dence detection of all three photon in the three different 
output modes. We will show that for a whole range of 
parameters of the 2D waveguide structure our system can 
suppress this coincidence event which corresponds to the 
original Hong-Ou-Mandel effect mm extended to three 
interfering photons. 


MODEL AND TIME EVOLUTION 

The investigated system consists of a 2D 3x3 waveg¬ 
uide array (three input modes and three output modes) 
with continuous evanescent coupling between the inner 
mode and the outer modes (see Fig. [l]). We assume 
identical single mode waveguides with uniform coupling 
throughout the waveguide array. We also assume that 


Hint = h (giaial + g 2 a 2 a\ + gla 2 a\ + g^a^ ■ (1) 

Each part of Eq. ([I]) stands for the annihilation of a 
photon in a certain mode and the creation of a photon 
in a neighbouring mode associated with the correspond¬ 
ing coupling strength g\/ 2 . Although g 1 and g 2 are real 
for waveguide systems we keep these parameters com¬ 
plex, since SU(3) Hamiltonians occur for many physical 
systems - for example N identical three level atoms in¬ 
teracting with external fields of different phases. We fur¬ 
ther add that the Hamiltonian that we use adequately 
describes the investigated waveguide system although it 
is not the most general SU(3) Hamiltonian. The calcu¬ 
lations can be done for the most general case, however, 
we will see that the Hamiltonian in Eq. 0 already leads 
to a number of very interesting interference effects by 
keeping the analysis more simple. 

In order to analyse the time dependent evolution of the 
system we switch to the Heisenberg picture. To simplify 

the calculation we define a vector a = ^dj, dj, ajj so 
that the time evolution is governed by 

d (° * °\ 

— a(t) = 1 ctf 0 g 2 \ a(t), (2) 

d V° 9*2 0/ 

V X 

M 

where the interaction time t is determined by the length 
of the waveguide. The equation can easily be solved us¬ 
ing the exponential ansatz V = which allows to 

rewrite the creation operators a, j (0) at time t = 0 in terms 

of creation operators d] (t) at time t via a(Q) = V ■ a(t). 
The solution of this equation yields the explicit form of 
the matrix V 


cos 2 (9) cos (G) + sin 2 (0) 

— icos(0) sin(G) e 1 ^ 

cos(0) sin(0)e i ^+^) (cos(G) - 1)\ 


-icos(0) sin(G)e“ 

cos(G) 

— i sin(0) sin(G)e ‘ v 

(3) 

cos(0) sin(0)e _ (cos (G) — 1) 

— isin(0) sin(G)e“ i¥3 

sin 2 (0) cos(G) + cos 2 (0) / 



where g\ -t = G cos (9) e 1 ^ and g 2 -t = G sin(0) e 1 v . From by V is periodic with respect to G and 9. The periodicity 
this expression it is easy to see that the evolution given in 6 arises due to the change to polar coordinates while 
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FIG. 2. Three examples for different indistinguishable quan¬ 
tum paths leading to the same output state 1300) in (a), 1021) 
in (b) and 1111) in (c) 

the periodicity in G is linked to the interpretation of G as 
overall scaling factor, i.e., for example to variations of the 
interaction lengths of the waveguide array. Specifically, 

for G = + | < 72 1 2 • t = 2tt ■ m, with m £ N, we 

get sin(G) = 0 and cos(G) = 1, so that V reduces to 
the identity matrix. This means that the system returns 
to its initial state, what is also known as self-imaging or 
revivals of multi mode interferences [ill 132] . Note that, 
as expected, the transformation V is unitary - this is 
because we are considering a lossless device. 

THREE SINGLE PHOTON INPUT 

Next we focus on an input state where at each input 
port a single photon is coupled into the waveguide at 
t = 0. The wave function of this state is given by \ipin) = 
a{(0)a2(0)dj(0) |0). Via the transformation matrix V we 
can easily calculate the general form of the output state 

3 3 3 

IV'out) = Y ■ Y V 2 m°L{t) ■ Y V 3nai(t) |0) , 

1=1 m =1 n =1 

(4) 

where V mn is the matrix element in the mth row and the 
nth column of V in Eq. As the transformation matrix 
V only acts on the creation operators of the three differ¬ 
ent modes in Eq. 0 a vacuum state is not transformed 
by the time evolution governed by the Hamiltonian of 
Eq. §. The general output state is a superposition of 
all possible distributions of the three photons among the 
three output modes with respective coefficients depend¬ 
ing on the explicit form of V. The general form of the 
output state reads 


where corresponding coefficients can be calculated by ex¬ 
plicitly expanding Eq. 0 or alternatively using a for¬ 
malism for linear optical networks involving permanents 
02 - 

A permanent of a matrix is equal to its determinant 
but without the sign of the permutation taken into ac¬ 
count. For a N x N matrix A it is given by 

N 

Perm A = £n (6) 

a j=l 

where the sum runs over all possible permutations er of 
the set {1,2,... ,N}. 

The coefficients Cu m of Eq. 0 can be expressed by 
permanents of matrices I/HT" 1 !, where k, l and m are 
the number of photons in the three output modes so that 
k + l + m = 3. Hereby V^ k,l,rn ^ is a 3 x 3 matrix and 
is constructed via the transformation matrix of Eq. 0 . 
It consists of k copies of the first column of V, l copies 
of the second column of V and m copies of the third 
column of V [43]. Dividing the permanent of by 

a normalisation factor yields the final expression for the 
coefficients 

PermV{ fe 4,m} 

° klm ~ ■ { ) 

One can show that the absolute value of these coefficients 
depends only on G and 9 but not on the phases 0 and ip 
of the coupling coefficients <q/ 2 which will only have an 
impact on the phases of the coefficients Ckim- 

Note that Eq. 0 only holds true for the input state 
11,1,1). However, as shown in [43], the formalism in¬ 
volving permanents can be expanded to arbitrary ini¬ 
tial states by additional consideration of the rows of the 
transformation matrix corresponding to the input state. 

The permanents can be understood as produced by 
the coherent superposition of indistinguishable quantum 
paths (excluding the normalization factor) leading to the 
same output state. To illustrate this we consider in the 
following three examples, starting with the coefficient 
C 3 oo- To calculate this coefficient, associated with the 
state 1300), one first has to construct the matrix I/f 300 ! 
containing three copies of the first row of the transfor¬ 
mation matrix V. It therefore reads 

(V n V u V-n\ 

U* 3 ’ 0 ’ 0 * = V 21 V 21 V 21 , (8) 

W V 31 V 31 ) 

where V nm are the corresponding matrix elements of V 
in Eq. 0 determined by the transition amplitude for 
a photon initially in mode n to exit in mode m. From 
Eq. 0 the permanent of I/f 30 °f can be calculated yield¬ 
ing 


IVw) —C 300 |3> 0,0) + CQ 30 10,3,0) + . • • 
+C 210 12,1,0) + ... + cm 11,1,1) 


(5) 


Perm V {3 ’°’ 0} = 6V u V 2l V :u . 


(9) 
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FIG. 3. Plot of the three photon coincidence probability 
|cm | 2 and the HOM contour, where the condition |dn| 2 =0 
is fulfilled. 


As a last example we calculate the coefficient cm cor¬ 
responding to a three mode coincidence event. As the 
matrix I/l 111 ! consists of a copy of each row it is equal 
to V and its permanent is calculated to be 

Perm =Vnl4 22 V33 + V12V23V31 + V13V21V32 

+V11V23V32 + V12V21V33 + V13V22V31. 

.. (12) 

It consists of six different terms each corresponding to 
a different indistinguishable quantum path leading to 
the same final state 11,1,1). For example, the first 
term corresponds to the case where all three photons 
exit the waveguide in the same mode they came in, 
the second term corresponds to the case where the pho¬ 
ton in the first/second/third mode switches to the sec¬ 
ond/third/first mode etc., as depicted in Fig. m^c). As in 
the second example the multiple indistinguishable quan¬ 
tum paths can interfere with each other. In particular, 
they can interfere in a completely destructive way. This 
configuration will be analyzed in the next section. 


This expression corresponds to the only possible quan¬ 
tum path where a photon in the first mode stays in the 
first mode and a photon of the second mode and a photon 
of the third mode also exits the waveguide array in the 
first mode as illustrated in Fig. da). As only a single 
quantum path appears for this coefficient no interference 
is observed in this case. 

The calculation of any other coefficient like, e.g., C 021 
for the state |0, 2,1), follows the same structure: First 
one constructs the matrix I/f 0 ’ 2 ’ 1 ! by taking two copies 
of the second row and one copy of the third row of V 


(V 12 V 12 
F {0 ’ 2 ’ 1} = F 22 H 22 

V 32 ^32 


V 13 \ 

V 23 ( 10 ) 

V 33 / 


and then calculates its permanent, yielding in this case 

PermF^ 0,2,1 ^ = 2(V’i 2 V r 2 2V'33 + V 12 V 23 V 32 . 

+ 1 ^ 13 ^ 221 / 32 ). 

Now three different indistinguishable quantum paths ap¬ 
pear which have to be added coherently as shown in Fig. 
§b)> what leads to interference effects. The three quan¬ 
tum paths correspond to the three possibilities how three 
photons entering the waveguide array initially in three 
different modes exit the array in the output state 1021). 


THREE PHOTON HONG-OU-MANDEL 
INTERFERENCE 

In the following we investigate the particular situation 
displaying a three photon Hong-Ou-Mandel (HOM) in¬ 
terference. In analogy to the original two photon Hong- 
Ou-Mandel experiment 111 m where two photons are 
never detected simultaneously at the two different out¬ 
put modes of a 50/50 beam splitter, the probability for all 
three photons leaving the waveguide at the three output 
ports vanishes if 


c 111 — 0. 


(13) 


To analyse the conditions for the three photon HOM 
interference, we have to calculate Cm explicitly. From 
Eqs. Q and (12) we find 


Cm =Vl ] V22 ^33 + V12V23V31 + V13V21F32 
+V11V23V32 + V12V21V33 + V13V22V31. 


(14) 


By inserting the expression for the various matrix ele¬ 
ments Vmn of Eq. ([3]) and solving Eq. (13) we can find 
an analytical expression for the HOM contour in the vari¬ 
able space (G, 9), where all states |Vw(G, #)) have a 
vanishing cm coefficient: 


9(G) = mr ± arcsec 



8 ± 


\/2 esc 4 ( / ) ^/sin 4 ( ^ ) (20 cos (G) +3(8 cos(2G)+4 cos(3G)+3 cos(4G)H-5)) 
3 cos(G)+2 


(15) 


As can be seen from Eq. (151 completely destructive range of the parameters g\ and gi. Note that for some 


three photon HOM interference can take place for a large 
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values of G these equations would result in a complex 
valued 6 and are therefore not considered as a solution. 
Fig. 0 shows a plot of the HOM contour, displaying in 
particular the periodicity in G and 0 discussed in Sect. 2. 

INTERESTING STATES ON THE 
HONG-OU-MANDEL CONTOUR 

Finally we investigate the states determined by the 
HOM contour. In the original HOM experiment a maxi¬ 
mally entangled state of the form oc |2, 0) — |0,2) is pro¬ 
duced at the output mm- Similar states can be found in 
the case of a three photon interference. Additionally to 
the condition cm = 0 we find that at certain points some 
coefficients Ckim of Eq. (0) will vanish as well, so that fur¬ 
ther terms are suppressed. Other coefficients will have 
the same absolute value, so that the states can be writ¬ 
ten in a compact form. We found three different kinds 
of states fulfilling this condition, which display entangle¬ 
ment between two and possibly three output modes. 

In Fig. ga) , where the HOM contour is shown, some 
coordinates are marked by points where one can find 
maximally bipartite entangled states. A closer investi¬ 
gation yields that for all these coordinates we find states 
of the form 


in Fig. 0(c). At the coordinates of the green dots (blue 
diamonds) it is evident that the outer mode mode l = 3 
(l = 1) is physically decoupled since the coupling param¬ 
eter g 2 {g i) vanishes. However, at the red crosses both 
coupling constants g± and g 2 are physically present but 
the inner mode l = 2 is effectively decoupled from the 
system and merely acts as a mediator between the outer 
modes. Therefore, it is more appropriate to speak of 
a two photon interference at these coordinates although 
three photons are present. The physical decoupling of 
modes has also been discussed in Campos [27] in the con¬ 
text of beam splitter devices. Note that all the discussed 
states are created in a deterministic way so that no post 
selection is required. 

In Figs. 0b) and [0] c) the coordinates of possible tri¬ 
partite entangled states along the HOM contour are dis¬ 
played. Note that we have three modes and each mode 
has four possible states corresponding to the occupation 
of 0, 1, 2 and 3 photons. We are thus dealing with higher 
dimensional entanglement of three qudits (d = 4). A 
complete classification of the classes of entangled states 
for qudits (d = 4) does not exist. However the structure 
of tripartite states generated for the case of Figs. 0b) 
and 0c) suggests three qudit (d = 4) entanglement. The 
form of the states of Fig. 0b) and their coordinates read 




IVW)_ ^2 ( 

|2„0 fc ) + |0„2 fe )) |1,) 

(16) 

l^out) = ^ ( 

J3j, 

o fc ) 

+ 

|0„ 

3*>) 

|0|) 






+ *( 

'\2j, 

lfc) 

+ 

1 lj, 

2 fc >) 

|0,) 

(17) 

with 

3 = 

1, k = 3 and 

l = 2 at the red crosses, 

with 

4 ' 








3 = 

1, k 

= 2 and l = 

3 at the green dots and 

with 

+ b( 

'\lj, 

Ofc) 

+ 

|0j, 

4)) 

|2,) 


3 = 

2, k 

= 3 and l = 

1 at the blue diamonds. 

Note 

2 ' 









that we neglected the global and relative phase factors 
of each state, as we just want to focus on the structure 
of the states (see the appendix for the exact analytical 
expressions for each coefficient and the corresponding co¬ 
ordinates) . All the three states have a similar form where 
one mode containing one photon is separable while the 
remaining two modes are in a maximally entangled state. 
Depending on the phase ip/^P °f the coupling coefficients 
gi /2 the relative phase between the non-separable states 
can be varied. A closer look at the transformation ma¬ 
trix V of Eq. (0) reveals that in these cases the abso¬ 
lute values of the matrix elements are given by |Vjj| = 1 , 
\V 3 k\ = \v kj \ = I Vjj\ = IVfcfcl = while all remain¬ 
ing matrix elements vanish. This means that the mode 
l is decoupled from the system as the photon initially 
entering the waveguide array in this mode will always 
exit the waveguide array in the same mode. However, 
the modes k and j are mixed with a 50/50 ratio corre¬ 
sponding to the original two photon HOM interference 


effect. This can also be seen from Eq. (14) where under 


the same conditions only two of the six quantum paths 
survive and interfere completely destructively. For 1 = 2 , 
the corresponding quantum paths are the two paths, one 
at the extreme left and the other at the extreme right, 


with j = 2, k = 3 and l = 1 at the red crosses, with j = 1, 
k = 3 and l = 2 at the green dots and with j = 1, k = 2 
and l = 3 at the blue diamonds. The form of the wave 
function of the states of Fig. 0(c) and their coordinates 
are given by 

IV’out) = /- ^2 13^, Ofc, 0/) + 10j, 3fc, 0;) + |0j, 0 fc , 3;)) 


3V2 
1 

76 
1 


+ ^(|l i ,0,) + |0 j ,l,>) |2 fc ) 


+ 


76 


(|lj,0/) + |0j,I*;}) |2,) 


(18) 


with j = 1, k = 2 and l = 3 at the red crosses, with 
j = 2, k = 1 and l = 3 at the green dots and with 
j = 3, k = 1 and l = 2 at the blue diamonds. As be¬ 
fore we neglected the global and relative phase factors of 
each state for simplicity (see the appendix for the exact 
analytical expressions for each coefficient and the corre¬ 
sponding coordinates). We have written the states in a 
way that a certain mode is always factored out in each 
term so that the entanglement between the two remain¬ 
ing modes is clearly visible; this suggests that the states 
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(a) (b) (c) 



FIG. 4. Plots of the HOM contour and the coordinates of the investigated states (red crosses, green dots and blue diamonds), 
which display bipartite entanglement (a) or possibly tripartite entanglement (b)/(c). 


are not separable and are good candidates for tripartite 
entanglement. 


CONCLUSION 

In conclusion, we investigated the dynamics of a 2D 
3x3 waveguide where the outer modes are coupled to 
the inner mode by evanescent coupling but not to each 
other. Beginning with three indistinguishable single pho¬ 
tons at the three input ports we showed that for a wide 
range of waveguide parameters this leads to completely 
destructive three photon interference, i.e. for these pa¬ 
rameters the photons will never leave the waveguide in 
three separate ports. This is a generalisation of the well 
know Hong-Ou-Mandel effect from two to three photons. 
Additionally, the produced output states consisting of 
three qudits (d = 4) exhibit highly interesting structures 
displaying bipartite or possibly tripartite entanglement. 
Clearly the entanglement of qudits is a topic for future 
studies. 
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For certain waveguide parameters G and 8 we found 
three interesting sets of states on the HOM contour for 
which the coincident event is suppressed and therefore 
Cm vanishes. However, only the general form of these 
states was discussed but not the explicit values of the 
coefficients. Here we present their analytical values: Ta¬ 
bles [I] to [V] contain the exact analytical expressions for 
the coefficients of Eq. (16) - ( |T8| ). Note that in Eqs. 
(Il7|) and (18) some signs depend on the value of nmod4. 


Therefore in the second part of Table |H|and in Table |TV| 
n is replaced by 4h + 0, 4fi + 1, 4n + 2 or 4h + 3, in which 
case all possible cases are considered. 
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TABLE III. Table of coefficients for equation (181. Coefficients, which are not listed, are equal to 0. 
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TABLE IV. Table of coefficients for equation (181. Coefficients, which are not listed, are equal to 0. 
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Table of coefficients for equation (|18|) . Coefficients, which are not listed, are equal to 0. 





























































